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Abstract 

We wish to study an application of Stueckelberg's relativistic quantum the- 
ory in the framework of general relativity. We study the form of the wave 
equation of a massive body in the presence of a Schwarzschild gravitational 
field. We treat the mathematical behavior of the wavefunction also around 
and beyond the horizon (r = 2M). Classically, within the horizon, the 
time component of the metric becomes spacelike and distance from the ori- 
gin singularity becomes timelike, suggesting an inevitable propagation of all 
matter within the horizon to a total collapse at r = 0. However, the quan- 
tum description of the wave function provides a different understanding of 
the behavior of matter within the horizon. We find that a test particle can 
almost never be found at the origin and is more probable to be found at 
the horizon. Matter outside the horizon has a very small wave length and 
therefore interference effects can be found only on a very small atomic scale. 
However, within the horizon, matter becomes totally "tachionic" and is po- 
tentially "spread" over all space. Small location uncertainties on the atomic 
scale become large around the horizon, and different mass components of the 
wave function can therefore interfere on a stellar scale. This interference phe- 
nomenon, where the probability of finding matter decreases as a function of 
the distance from the horizon, appears as an effective gravitational repulsion. 



1 Introduction 



In Newton's classical mechanics and in quantum mechanics, one makes use 
of a global time that has causal meaning. In standard non-relativistic quan- 
tum mechanics, time is interpreted as a causal parameter, where, for each 
value of the parameter, the quantum states are coherent. The manifestly 
covariant quantum Stueckelberg formalism is based on the idea that there is 
an invariant parameter r of evolution of the system; wave functions, as co- 
variant functions of space and the Einstein time t form a Hilbert space (over 
R A ) for each value of r. Thus, there are two types of time, one transforming 
covariantly, and the second, a parameter of evolution pQ. 

J. A. Wheeler was one of the first physicists to understand the physical 
difference between both kinds of time. Together with his student, Richard P. 
Feynman, they showed the significance of both the retarded and the advanced 
parts of the particle propagator for the correct description of direct particle 
interaction [2] (in the 5D electrodynamics associated with the Stueckelberg- 
Schrodinger equation, the r-retarded propagator contains both t-retarded 
and advanced components [3]). The symmetry between past and future in 
the prescription of the fields as a consequence of the QED theory, leads to 
the understanding that a causal parameter in our experience, i.e., an in- 
variant universal time, may not be totally correlated with the Minkowski or 
Einstein time, in which interactions notice past and future simultaneously. 
This means that the system admits quantum superposition of states in the 
total 4D picture where time acts as another spatial dimension and events can 
move forward and backward in time t. Later, with his famous Delayed Choice 
thought experiment [4], Wheeler showed that a measurement may intervene 
in the static 4D picture and alter the whole 4D picture in its entirety. The 
intervening measurement not only changes the "present" state (an incident 
that would disturb unitary evolution and other conservation laws), but actu- 
ally transforms the entire "history" of the 4D state and therefore maintains 
unitarity; four-momentum is conserved [5]. 

In relativity, the time of an event as measured in the laboratory is subject 
to variation according to the velocity of the apparatus related to the trans- 
mitting system and may as well be affected by forces (such as gravity). A 
relativistic quantum theory should therefore incorporate time in a manifestly 
covariant manner and also permit definition of a global causal parameter to 
generate evolution of the 4D states. 

To describe the dynamical evolution of such a system, Stueckelberg and 
Horwitz-Piron [61 [7] introduced an invariant parameter r, a "world time", 
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coinciding with the Newtonian world time, accounting for classical as well 
as quantum relativistic evolution. The world time r provides a parameter 
that labels the dynamical evolution of the covariant system. For free on shell 
motion of a single particle, the Einstein proper time can be taken equal to 
the world time. 

We shall make use of this Manifestly Covariant Quantum Theory (MCQT) 
introduced by Stueckelberg in 1941 and generalized by Horwitz and Piron 
[7] in 1973. The theory has the structure of Hamilton dynamics with the 
Euclidean three dimensional space replaced by four- dimensional Minkowski 
space. Since all four components of energy-momentum are kinematically 
independent, the theory is intrinsically "off-shell". 

The self-interaction problem of the relativistic charged particle has re- 
cently been studied, where it was shown that the radiation field is associated 
with excursions from the mass shell [H [101 El HH 02] • In the limit that the 
motion maintains a state very close to the mass shell limit, the equation re- 
duces [5] to that of Dirac [13] (the Abraham-Dirac-Lorentz equation). The 
recent experiment of Lindner et al [H] , showed the phenomenon of interfer- 
ence in time, an essential property of quantum mechanical wave functions 
coherent in time, as in the Stueckelberg theory. The results were predicted 
by Horwitz and Rabin in 1976 [15], and the significance of this important 
experiment was discussed in [TO] . 

Another result which concerns the Stueckelberg MCQT formalism which 
was pointed out recently [17] , is that by using the metric tensor in the kine- 
matical terms of the Stueckelberg-Schrodinger equation one can obtain clas- 
sical general relativity in eikonal approximation. Since the eikonal approxi- 
mation lowers the dimension of the differential equations describing the fields 
by one [IS], the eikonal approximation to the five dimensional Stueckelberg 
quantum equation in a curved spacetime, characterized by a metric tensor 
g^v, results in the four dimensional Einstein geodesic equations. 

When dealing with a gravitational potential represented by the metric 
tensor within the covariant theory, the theory remains on-shell and can be 
shown to have a form similar to the Klein-Gordon (KG) equation on a curved 
spacetime for each value of the mass (on a continuum). 

In this paper we study an application of Stueckelberg's MCQT in gen- 
eral relativity. We deal with a simple case and compare the results to those 
expected in classical general relativity. We study the form of the wave equa- 
tion of a test particle in the presence of a Schwarzschild gravitational field, 
assuming that the source is massive enough to ignore changes in the metric 
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caused by the smaller mass of the test particle. J. Makela [19] has consid- 
ered the black hole as an atom in the framework of quantum theory, but did 
not discuss the spatial distribution of the wave functions. Since the Stueck- 
elberg formalism describes the evolution of the wave function of the body 
according to the invariant world time of evolution, we can treat the math- 
ematical behavior of the wavefunction also around and beyond the horizon 
(r = 2M). Within the horizon, the interval classically effectively changes 
its signature and becomes spacelike; the distance from the origin singularity 
becomes timelike, but the description according to world time enables us to 
study the behavior of the evolution of the wavefunction in an an absolute 
sense. 

We find that within the black hole horizon, the expectation value of the 
distance from the R = singularity has a strong gradient towards the hori- 
zon. This result can only be explained quantum mechanically, since classi- 
cally, the particle should move towards the origin. Interference effects ev- 
idently induce results which are very different from the classical behavior. 
The interference makes a difference in the results only when the particle has 
a wavelength of stellar scale, A ~ M. In the neighborhood of the horizon, the 
wavefunction is evidently spread, in analogy to the action of tidal forces. This 
result means that for certain wave-packets, there is "gravitational repulsion" 
that prevents the test particle from falling towards the R = singularity, 
and maintains it near the horizon in a manner that partly depends on the 
particle's angular momentum. Therefore, this quantum gravitational model 
predicts that the test particle will move to the r = 2M shell of the black 
hole, which means that matter in the black hole should accumulate on the 
interior of the horizon. As we have indicated above, in the eikonal (ray) 
approximation, the quantum solutions flow as in the classical limit, and our 
interpretations are therefore a valid description of the expected observed be- 
havior of a particle as a result of the application of the quantum theory. 

The dynamics associated with the matter repulsion within the black hole, 
causes the Einstein time to run in the opposite direction of the world time. 
This effect is well-known in the Stueckelberg-Feynman diagram of electron- 
positron annihilation, where the positron is interpreted as an electron going 
"backwards" in time. This effect around the horizon will cause matter within 
the black-hole to have properties of antimatter. This model may explain some 
part of the matter-anti-matter asymmetry observed in the universe, if there 
are sufficient black hole surfaces. Using covariant quantum theory, therefore, 
may help provide a new understanding of black hole physics and cosmology. 

This paper is in memory of John Archibald Wheeler, who passed away 
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the week this paper had started being written. In this paper we combine 
some of his ideas regarding time, together with the physics of the "black 
hole" , a term which he came up with and brought to worldwide attention. 



2 The Eikonal Approximation 

We start with the Stueckelberg-Schrodinger (Horwitz-Oron) equation [201 
EE]: 

where m is an intrinsic property of the particle with dimension mass. 

The Stueckelberg Hamiltonian is written with a non-trivial metric type 
function g^ u on spacetime in the quadratic kinetic term. 

We will show hereafter, that for the eikonal (semi-classical) approxima- 
tion, these equations lead to Einstein's geodesic flow on a curved manifold; 
general relativity then appears as an emergent phenomenon [20J. 



In the eikonal approximation for the 5D equation above, we assume a \I/ 
such that: 

V(x, t) = A(x)e- i( ^ T -^ s{x)) (2) 

which gives: 
2m 2m m l 

(3) 

In the eikonal approximation, it is assumed that K (k = p^ = — m 2 , the 
dynamical measured mass squared; we use the signature (-, +,+,+) in the 
local flat space) is large compared to the square of the second derivative of 
S; the dynamical evolution of the system in r is effectively frozen, and the 
theory reduces to a four dimensional eikonal form: 

d»Sd v S 

g»u — = -l 4 

m z 

In the optical analogy of the eikonal approximation (22], the functions 
S(x) are the Fresnel surfaces of rays, and d^S = is the momentum in 
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the direction of the propagation of the phase surface. If we identify the 
momentum p^ as mx M where the dot corresponds to differentiation in r, this 
equation becomes: 

— g iiv dx p 'dx v = dr 2 (5) 
the equation for the invariant line element of Einstein. 
The Schrodinger current associated with (TjQ) is [2"2"j [23]: 

Jr(x) u = JL(** gitt ,di>* T - ^ T g^K) (6) 

The Fresnel surface of the system's dynamics, analogous to the optical 

case is |22|: 



K = g^V u + m 2 = (7) 

It is clear that dK/dp^ is in the direction of the eikonal form of j£. 
This implies that K is the evolution operator for the dynamical flow of the 
particles which correspond to the Fresnel rays. K is therefore, the covariant 
Hamiltonian of the system. It then follows from the Hamilton equations that 
the flow is geodesic, where g^ u is the metric [2"5] . 

We note that according to the Stueckelberg theory, ds/dr = m/m, where 
^ 2 = 9iiuP fl p u ' i and this change of variable in §5§ brings the line element to 
the form using proper time ds with a factor ih/m, said to be "on shell" if 
this quantity is unity. 



3 The Stueckelberg Formalism in a Schwarzschild 
gravitational field 



For the Schwarzschild metric we have (we consider the purely radial case for 
which = 0, and take G=l): 



( 


1 





^ 2M 




r 




2M 




r 










\ 












o r - 



\ 



esc 2 8 



and also: 



'—g = r 2 sin 9 



(9) 
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After substituting the metric into equation <^ and separating variables, we 
find that the solution to the problem is of the form: 

M/ T ( x ») = e-^-^-^P^R^r) (10) 

where k = 2mk has the dimensions of mass 2 , and equals m 2 on the particle's 
mass shell. 

Substituting this form into equation ([I]), we get an equation for R K}U)j i(r): 



R>)+ ^-M) f^W ^ + /(/ + l) \ 

i j + r(r-2M) { ) + \(r-2M) 2 r(r-2M) ) K ) K ' 

and for P^ a {9): 

P"{6) + cot 0P'(6) + + 1) - -^-) P{6) = (12) 

\ sin V J 

3.1 The form of the gravitational potential in the weak 
gravity limit (r >> 2M) 

For the far gravitational field, where we take r » 2M and therefore neglect 
high orders of — , after substituting k — > m 2 , the radial equation has the 
form: 



„.,. , 2„. /N . „ „ 1(1 + 1) 2M(2uj 2 -m 2 )^ r . . , 

R"{r) + -R'{r) + {u 2 - m 2 - y ' + ^ L )R{r) = 13 

If we divide the whole equation by —2m we find that this equation has 
the exact form of the Schrodinger central potential equation, where the term 
represents the repelling centrifugal potential. The Centrifugal term shall 
be omitted hereafter, since we shall keep only the first order terms of -. 

Since, relativistically, u is the energy of the particle including its mass 
and since we are not interested here in tachyonic solutions, we take u 2 —m 2 to 
be positive. Therefore, as opposed to the nonrelativistic hydrogen problem, 
the eigenvalues of interest are positive. The tachyonic solutions exponentially 
decrease at infinity and therefore don't add any relevant amplitude to the 
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far field solutions; however, one could argue that they should be taken into 
consideration near the horizon. 

The effective central potential for this equation has the form: 



M(2uj 2 -m 2 ) M . 2p 2 N mM . 9n . . 

V(r) = m+— )« (l + 2v 2 ) 14 

mr r m r 

The effective gravitational mass is altered by a small factor of 2v 2 . Moving 
back to ordinary units, the effective gravitational mass is: \x = m(l + ^-). 

The solution to the radial equation is of the form: 

R{r) = clJi (2yggg + C2Y, (2^M¥) (ig) 



where J n (z) and Y n (z), give the Bessel function of the first and second 
kind accordingly. 

The fact that the energy eigenvalues of the equation are positive means 
that the wavefunctions, as single momentum modes, are not square inte- 
grable. We consider therefore wave packets which include many momentum 
modes. For the far fields, our massive object may be localized very tightly, 
and may have an uncertainty in its location which is very small compared 
to the scale of the gravitational distances. This situation, where the object's 
wavelength doesn't play a role in the physics of the problem, changes when 
the metric causes the wavelength (or equivalently the location uncertainty) 
of the object to be of relevant scale, which will cause interference phenomena 
to be crucial. 

We now discuss predictions close to and within the black hole horizon, and 
in a following section show the results of a complete and exact computation 
of the wavefunction in the entire interior region. 



3.2 Solutions for the radial part in special cases 

Taking: 

B(r) 

R{r) = r A\m (16) 

yfr^/r — 2M 

and substituting in equation ( ITT]) , we get for B m ^,i{r) a Schrodinger-like 
equation of the form: 
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_ r(r - 2M) gg + 1) + mV) - M 2 - 
1 ; (r-2M)V 1 ; 1 J 



3.2.1 Solving the radial equation when approaching the horizon 

A general solution of the radial equation seems to be very difficult to achieve. 
Therefore, we solve the equation for those regions where simplifying assump- 
tions permit an easier calculation. We know that such assumptions contribute 
to a certain amount of inaccuracy, and therefore we can treat the results only 
qualitatively and not quantitatively. The solutions give us a good idea of the 
direction in which the body will move. When the body moves to other re- 
gions, we must take into account the other parts of the potential, and the 
calculations become more complicated. 

We first study the behavior of particles near the horizon (r — > 2M + ). 
Much work has been done on the behavior of the field equations when ap- 
proaching the Schwarzschild horizon from above, for example, the most gen- 
eral solution of the radial equation in an asymptotic form, under the assump- 
tion of independency of time (which is a good assumption at the horizon), is 
given in [24]. 

Expanding the potential around r — > 2M + , the potential in equation ([T 
takes the form: 



to 



2 



2u 2 - m 2 1(1 + 1) + ~ 



u w = -# — + ^kr < 18 > 

where d = r -^f 0. 

This is clearly a very strong attracting potential (apart from the relatively 
neglected weak repelling centrifugal part). 

Solving the radial equation in this region we obtain for the dominant 
parts: 

R{d) = a x d~ 2iMw + a 2 d 2iMw (19) 

Finding the expectation value of r — > 2M + , with normalization, we need 
to compute: 
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f™ +X 4irr 3 R(r)R(r)*dr 
l2M +X 47rr 2 i?(r)i?(r)* dr 



< r >= n „\ [1 { , (20) 



where A is the uncertainty in the location of the particle on the R axis or, 
equivalently, the width of the wave packet associated with the wavelength of 
the particle. 

Changing variables around the horizon we get: 



M (aia 2 ((3e - 4iMu)e 8 * Muj + AiMu) + 16M 2 u 2 e 4iMuJ ) 
<r>= aia2 (( e _ 2iMcu)e 8iMuj + 2iMu) + 8M 2 u 2 e 4iM " ^ 2M 

where e = ^ << 1- 

What can be seen from the result is the fact that the particle is very 
strongly captured by the horizon. The uncertainty in its location is now 
related to the phase in the term e tMuJ which becomes exponentially small and 
concentrated near the horizon through normalization of the wavefunction. 



3.2.2 Solving the radial equation at r — > 2M 

Taking the same potential from equation (Tl8l) and solving the radial equation 
for r — > 2M~ we obtain (only the first term of (1191) applies to the interior 
solution): 

R(d) = ai d~ 2iMuJ (21) 

where d = _> . 

Finding the expectation value of r in a small neighborhood of the horizon, 
i.e., at r = 2M — p, where p « 2M we need to compute: 



< r >= f™:£4irr>R(r)R{rydr 
J 2 \ A !: p P ^^r 2 R(r)R(rydr 

Changing variables around the horizon we get: 

2A 2 

<r>=2M~p + 
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It can be seen that within a small distance from the horizon at the point 
r = 2M — p, the expectation value tends towards r = 2M. This result means 
that effectively, a particle that is at the horizon, will most probably stay 
there and not fall into the center. 

Using A in the expectation value calculation indicates a width for com- 
puting the local probability. Since we are using A as the effective spread of 
the particle's wave function (and beyond that we have a zero probability of 
finding the particle) the calculation is a good approximation. Of course, if 
A is spread beyond the validity of our approximation, this assumption isn't 
justified. 

The reason that the interacting particle behaves differently quantum- 
mechanically than is expected classically, is the fact that the expectation 
value takes interference into account. A similar phenomena is the reason for 
which the electron in its ground state doesn't fall into the atom. One can 
think of this as the result of the quantum effect of interference (see also [19]). 



3.2.3 Solving the radial equation when approaching r=0 

The equations are well defined within the horizon. Therefore, it is of interest 
to study the equation when r — > 0. 

When r — > (r << M), the potential in equation f|T7|) takes the form (a 
repulsive effective potential): 



U(r) - 2/(/ + 1) + 1 - — (23) 



We get for the radial equation the solution: 



R(r) = a 1 + a 2 log (^) (24) 

Finding the expectation value of r very close to the origin, at r = p, 
where p « 2M we need to compute: 



J^A7rr 2 R(r)R(r)*dr 



after some series expansion we get: 
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0.0 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ' A 
10 20 30 40 50 

Figure 1: The amount of " repulsion" from the particle's location (in A units) 
as function of the distance from r = (in A units) 



2A 2 p 

< r >Ri p H ■ — 

P A 2 + 3p 2 



The point we want to show in (1251) is, that the expectation value for 
a "particle" located initially at r = p is greater than p. Of course the 
expectation value is still small and doesn't "explode" at r = 0, and also the 
"repulsion" (i.e., the distance between the location p and the expectation 
value) becomes smaller and smaller as you draw away from r = (Figure 
[1]); however, this is a consistent result, and will cause matter eventually not 
to be found close to the origin even if it starts there. Therefore we can 
summarize, that, although the potential shows a very strong attractive force, 
the quantum mechanical equation implies a repulsive behavior. 

When p — > A, becomes smaller (if for instance the black hole was created 
with mass around the origin), the effect of repulsion becomes even more 
dominant. At the limit where p < A, the expectation value becomes: 

_ j P+X 4Tcr 3 R{r)R{r)*dr _ 3(A + p) 
~ f? +X 4nr 2 R(r)R(rydr~ 4 

The result we obtain shows that mass effectively moves from within the 
black hole to the horizon. 
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3.3 Exact computation of the wavefunction in the en- 
tire interior region 

In the previous section, we have computed approximations to the wavefunc- 
tion in three short ranges of radius, one near the origin, one near the horizon 
from within and a third near the horizon from without. However, in order 
to validate our conclusions from the previous section, we shall perform a 
calculation with a single and globally defined wavefunction that allows com- 
parison of the probability residing in separate regions. We shall, therefore, 
solve ffTTl) numerically for particular selected values of k and u and with 
proper boundary conditions, and give a full picture of the wave function in 
the whole (internal) region. 

3.3.1 Constructing the initial wavefunction 

For the numerical evaluation we shall choose: M = 100. Our particle's 
energy shall be described by a gaussian distribution function around a main 
energy of uq = 2.5 with a frequency width of cr w = 0.2. The mass of the 
particle shall be taken to be constant (with no distribution) over the whole 
wavefunction in order that the numerics fit the KG equation for a specific 
mass; m = \fn — ^ = 1.25. We shall also take the trivial case where / = 0, 
and so the wave function is independent of 9. 

Since t is spacelike within the interior, our initial conditions for the wave- 
function are a well localized gaussian in t (and in u) with width a = M/10 = 
10 (a w = 0.2 in the frequency domain). The gaussian is constructed numer- 
ically using 100 eigenfunctions with different frequencies around the central 
frequency ojq, with a width of A w = 1, such that the minimal frequency 
ojq — A w , is still greater than m and therefore we avoid tachionic modes in 
the solution. 

The wavefunction, as function of the spatial coordinates t and r, is there- 
fore described by the superposition of the states according to: 

$(r,t)= R "i(. r ) e ~ imt ( 26 ) 

i=1..100 

The first boundary condition on the evolution of $(r, t) is the gaussian 
as function of t, shown in Figure [2j when taking all the 100 frequency modes 
at r = M. 
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Figure 2: The boundary condition for the wavef unction at r = M 

The second boundary condition on the evolution of states shall be R' w . (ro = 
M) = 0. 

3.3.2 The "time-like" evolution of the wavepacket 

We shall now let the function evolve in V, which is thought as the timelike 
axis. 

After solving numerically equation (TTTjl . we can see in figure |3] the evolu- 
tion of the wavefunction $(r, t), in the interior region < r < M. 

Because of the symmetry of the boundary conditions around t = 0, the 
function separates into two wave functions, one propagating in the t direc- 
tion and the other propagating in the —t direction as r decreases. When 
approaching to r = 0, the propagation "freezes" in t and the only thing that 
changes in the evolution of the wavefunction is the amplitude which decreases 
as described in our analytic solution according to equation (j2lj) . The agree- 
ment between the numerical solution and our analytic solution near r = 0, 
is shown in figure HI 

The evolution of the wavefunction from the horizon to r = M is shown 
in figure [5j It can be seen that the boundary conditions cause two major 
parts of the wave function at the origin to interfere at r = M and create our 
gaussian. 

In figure [6j cuts according to r of the 3D probability density function are 
shown. 
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Figure 4: The agreement between the analytical to the numerical solution 
when approaching r = is shown. The purple dashed line is the analytic 
solution close to r = 0. 
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Figure 6: $(r, £)$*(r, t)r 2 - 2D cuts from the horizon to r = 
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The most surprising fact in our solution is the ongoing "disappearing" 
of the probability density function when approaching r = 0. According to 
the classical interpretation it looks as if the particle disappears in it's future, 
something that is of course impossible. The dependence of the wavefunction 
in t is therefore the only part in the wave function which is "non-localized" 
and therefore is the only variable that can be interpreted to be the 'time- 
evolution' of the particle. 

Since the numerical wavefunction we have created isn't localized in r, it 
is hard to normalize it around the horizon and therefore, it is impossible to 
perform an accurate expectation value for r. However, integrating over the 
probability density close to the horizon gives a bigger probability than the 
integration around r = M as opposed to the wave probability density func- 
tion at r=0 which is zero. The numerical solution therefore, also strengthens 
our previous analytic results that the particle is preferably found near the 
horizon. 



4 Discussion and Summary 



The quantum model we have used in this paper, based on the Stueckelberg- 
Schrodinger equation, suggests new physics regarding the prediction of a 
repulsive gravitational behavior within the black hole. This result can only be 
explained quantum mechanically. We have solved the r independent equation 
for r — > oo, the limit r — > 2M ± , and r — t- 0, and find that, as expected, the 
test particle at 2M + becomes trapped, but inside the horizon, it does not fall 
to the center, but (as for other initial conditions within the horizon) tends 
to be found around the horizon. Due to the similarity of the Stueckelberg 
theory to the structure of the KG (Klein-Gordon) equation, apparently, one 
could have used the KG equation and get the same results. As a matter of 
fact, equation (jTTj) is actually the radial part of the KG equation (if one, 
considers it as a wave equation). The KG equation, however, as pointed out 
by Newton and Wigner [22] does not provide a wave function from which 
one can make local conclusions about the distribution of matter, which is 
our purpose here. The Stueckelberg theory, however, which is intrinsically 
off-shell, has the property that it provides a quantum theory with the correct 
properties of locality [26] . It is for this reason that we studied this question 
in the framework of Stueckelberg. The results of our computations in this 
framework have a clear interpretation whereas from the point of view of KG 
theory, they are difficult to interpret. This is because Stueckelberg's "world 
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time" doesn't necessarily propagate with V within the black hole, while in 
the KG equation V is the time-like axis where the particle must propagate 
towards its "future" at r = 0, and must not "disappear" as is seems to do in 
a detailed analysis. 

We would like to thank Prof. Amos Ori from the Physics department of 
the Technion Institute for his valuable contributions to this article and espe- 
cially for his help in constructing a numerical analysis model that corresponds 
to a correct General Relativity approach. 



18 



References 

[1] L.P. Horwitz, R.I. Arshansky and A. Elitzur, Found. Phys. 18, 1159 
(1989). 

[2] J. A. Wheeler and R. P. Feynman, Rev. Mod. Physics 21, 425 (1949). 

[3] I. Aharonovich and L.P. Horwitz submitted; arXiv math-ph/0910.3343. 

[4] J. A. Wheeler, in Quantum Theory of Measurement (ed. Wheeler J. A. 
and Zurek W.H, Princeton U.P. Princeton), 182 (1983); J. A. Wheeler, 
in Mathematical Foundations of Quantum Theory (ed. Marlow A.R Aca- 
demic New York.), 9 (1978). 

[5] The Role of Momentum Transfer in Welcher-Weg Experiments, D. M. 
Ludwin and Y. Ben-Aryeh, Foundations of Physics Letters 14(6), 519- 
528 (2001). 

[6] E. C. G. Stueckelberg, Helv. Phys. Acta 14, 322, 588 (1941). See also: 
V.A. Fock, Phys. Z. Sowjetunion 12 (1937) 404; J. S. Schwinger, Phys. 
Rev. 82, 664 (1951); R. P. Feynman, Rev. Mod. Phys. 20, 367 (1948); 
R. P. Feynman, Phys. Rev. 80, 440 (1950) 

[7] Relativistic dynamics, L.P. Horwitz and C. Piron, Helv. Phys. Acta 
46, 316 (1973). See also, E.C.G. Stueckelberg, Helv. Phys. Acta 15, 23 
(1942). 

[8] O. Oron and L. P. Horwitz, Phys. Lett. A 280, 265 (2001). 

[9] N. Katz, L.P. Horwitz and O. Oron, Discrete Dynamics in Nature and 
Society (DDNS) 1, 179 (2004). 

[10] O.Oron and L.P. Horwitz, Phys. Lett. A 280, 215 (2001). 

[11] L.P. Horwitz and O. Oron, in "Quantum Gravity Research Trends, vol. 
250, Chap. 5, p. 143, ed. A. Reiner, Novs Press, Hauppage (2006). 

[12] A. Roitgrund and L.P. Horwitz, "Simulation of the Radiation Reaction 
Orbits of a Classical Relativistic Charged Particle with Generalized Off- 
Shell Lorentz Force," (to be submitted for publication). 

[13] PA.M Dirac, Proc. Roy. Soc. London, Ser. A 167, 148 (1938). 

[14] Lindner et al. Phys. Rev. Lett. 95, 040401 (2005) ;arXiv 
|quant-ph/0 503165 



19 



Relativistic Diffraction, L.P. Horwitz and Y. Rabin, Lettere al Nuovo 
Cimento 17, 501 (1976). 

On the significance of a recent experiment demonstrating quantum in- 
terference in time, Lawrence Horwitz, Phy. Lett. A, 3551-6 (2006). 

Eikonal approximation to 5D wave equations and the 4D spacetime 
metric, O. Oron and L.P. Horwitz, Found, of Phys. 33 (8) 1323-1338 
(2003)(hep-ph/020518). 

M. Visser, Class, and Quantum Grav. 18, 3595 (2001). 
J. Makela, Found, of Phys. 32, 1809 (2002). 

Relativistic Brownian motion in 3+1 dimensions, O. Oron and L.P. Hor- 



witz, sub. Jour. Math. Phys. Qmath-ph/0312003[ TAUP 2752-03) 



Eikonal Approximation to 5D Wave equations as geodesic motion in a 
curved 4D spacetime, O. Oron and L.P. Horwitz, Gen. Rel. and Grav. 
(TAUP 2751-03), Nova Science Pub. Hauppage, 109-122 (2004). 

M. Kline and I. W. Kay, Electromagnetic Theory And Geometrical Op- 
tics, John Wiley and Sons NY. Am. (1965), Ch. II,III 

C. W. Misner, K. S. Thorne and J. A. Wheeler, in Gravitation, W. H. 
Freeman, San Francisco (1973). 

S. M. Christensen and S. A. Fulling, Physical Review D 15, 2088 (1977). 

T.D. Newton and E.P Wigner, Rev. Mod. Phys. 21, 400 (1949). 

Relativistic dynamics, L.P. Horwitz and C. Piron, Helv. Phys. Acta 46, 
316 (1973). 



20 



